A computational study of cycloaddition reactions of d8 metal tetroxide (iron, ruthenium, osmium) complexes with C60.

The potential energy surfaces of the cycloaddition reactions MO(4)(NC(5)H(5))(2) + C(60)→ MO(4)(NC(5)H(5))(2)(C(60)) (M = Fe, Ru, and Os) have been studied at the B3LYP/LANL2DZ level of theory. It has been found that there should be two competing pathways in these reactions, which can be classified as a [6,5]-attack (path A) and a [6,6]-attack (path B). Our theoretical calculations indicate that, given the same reaction conditions, the cycloaddition reaction of C(60)via [6,6]-attack is more favorable than that via [6,5]-attack both kinetically and thermodynamically. This is in good agreement with the available experimental observations. A qualitative model, which is based on the theory of Pross and Shaik, has been used to develop an explanation for the barrier heights. As a result, our theoretical findings suggest that the singlet-triplet splitting ΔE(st) (= E(triplet)- E(singlet)) of the d(8) MO(4)(NC(5)H(5))(2) and C(60) species can be a guide to predict their reactivity towards cycloaddition. Our model results demonstrate that the reactivity of d(8) metal tetroxide cycloaddition to C(60) decreases in the order FeO(4)(NC(5)H(5))(2) > RuO(4)(NC(5)H(5))(2) > OsO(4)(NC(5)H(5))(2). In consequence, we show that both electronic and geometric effects play a decisive role in determining the energy barriers as well as the reaction enthalpy.